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SPECTRAL VALIDATION OF THE WHITHAM EQUATIONS FOR 
PERIODIC WAVES OF LATTICE DYNAMICAL SYSTEMS 

BUGRA KABIL AND L.MIGUEL RODRIGUES 


Abstract. In the present contribution we investigate some features of dynamical lattice sys¬ 
tems near periodic traveling waves. First, following the formal averaging method of Whitham, 
we derive modulation systems expected to drive at main order the time evolution of slowly mod¬ 
ulated wavetrains. Then, for waves whose period is commensurable to the lattice, we prove that 
the formally-derived first-order averaged system must be at least weakly hyperbolic if the back¬ 
ground waves are to be spectrally stable, and, when weak hyperbolicity is met, the characteristic 
velocities of the modulation system provide group velocities of the original system. Historically, 
for dynamical evolutions obeying partial differential equations, this has been proved, according 
to increasing level of algebraic complexity, first for systems of reaction-diffusion type, then for 
generic systems of balance laws, at last for Hamiltonian systems. Here, for their semi-discrete 
counterparts, we give at once simultaneous proofs for all these cases. Our main analytical tool 
is the discrete Bloch transform, a discrete analogous to the continuous Bloch transform. Nev¬ 
ertheless, we needed to overcome the absence of genuine space-translation invariance, a key 
ingredient of continuous analyses. 
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1. Introduction 

A common strategy to tackle the analysis of large time behavior of dynamical systems focuses 
on coherent structures, often playing the role of organizing centers for the time evolution. We 
follow here this line of investigation for lattice dynamical systems and restrict our attention to 
periodic traveling waves. 
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By definition, evolution of lattice dynamical systems is continuous in time but discrete in 
space. We consider here particular instances of those where the lattice is one-dimensional and 
thus may be assumed to be Z, and the evolution obeys a differential equation 

it) = T{U{t)) 

where the unknown U associates to any time t G R an R'^-valued sequence U{t) G (R*^)^ 
{d G N* being a fixed dimension) and is a vector-field on (R*^)^ that preserves the linear 
space of finitely-supported sequences and acts smoothly on it. Alternatively one may view the 
evolution as given by an infinite number of differential equations on R'^-valued functions coupled 
by a finite-range interaction. We expect that our results could actually deal with nonlocalized 
interaction satisfying certain short-range decaying assumptions but we choose to ignore those 
here to keep technicalities as low as possible. 

There are at least two ways in which such systems arise. They may come directly from 
modeling at a discrete level, as for neural networks, granular media, crystals, biological molecules, 
optical waveguides, chains of coupled oscillators... Alternatively they may emerge from the 
discretization in space — often called semi-discretization — of systems of partial differential 
equations, that generate fully continuous dynamical systems. In the latter case, mark that 
although a study of dynamical evolutions generated by fully discrete schemes including boundary 
conditions would be of much more direct practical use, such a study appears as a daunting task 
and it is a common belief that still relevant pieces of information are obtained from investigation 
of extended semi-discretized schemes. The reader looking for panoramas on lattice dynamical 
systems is referred to [?, ?, ?, ?, ?]. Specially designed technical tools of wide application range 
may be found in [?, ?, ?, ?, ?, ?, ?, ?, ?, ?, ?, ?, ?, ?, ?, ?]. Especially the latter references 
include key-ingredients to prove existence of traveling waves for various class of systems, a fact 
that is taken as an assumption in the present piece of work. 

We are interested here in a particular area of the dynamics, namely what occurs in neighbor¬ 
hoods of periodic traveling waves. Compared to other coherent structures of lattice dynamical 
systems — such as fronts, kinks, pulses, shocks, solitary waves or breathers —, periodic waves 
have received relatively less attention. In particular, while the stability analysis of some of the 
former patterns seems to have reached some maturity — see for instance [?, ?, ?, ?, ?, ?, ?, ?] 
—, the authors of the present paper are not aware of even a single example of a comprehensive 
stability study for periodic waves^ of some lattice dynamical system. Nevertheless, the interested 
reader may benefit from looking at [?, ?, ?, ?, ?, ?, ?, ?, ?, ?, ?, ?, ?, ?, ?, ?, ?, ?, ?], that are 
at least partly dedicated to periodic waves of lattice dynamical systems, and mostly focusing on 
proving their existence. 

To set our precise framework, let us specify that we say that, for systems as above, a traveling 
wave, sometimes called a uniformly sliding solution, — that is, a solution U with special form 
U(t) = {u{j — ct))j^z for some speed c G R and some profile tt : R ^ R*^ — is said to be periodic 
if the corresponding profile u is itself periodic. In this case we rescale phase variable of the profile 
to ensure that the profile period is one. Explicitly any periodic traveling wave U is written as 
U{t) = {u{k j + ujt))j,=z with a one-periodic profile u. This brings out spatial wavenumber 
k G R* and time frequency u = —kc. Mark that up to now we are not assuming that the 
spatial wavenumber k is rational. When its is not the case the spatial periodicity of the profile 
is not a priori easily observed on the solution itself, though when standing waves are discarded, 
that is when restricting to c 7 ^ 0 this periodicity is apparent in time, the solution being l/uj- 
periodic in time. We refer the reader to [?] for further elaboration on this comment including 
fascinating illustrating pictures. This already stresses that, in contrast with analogous issues 
for continuous dynamical systems, analysis of the behavior near standing waves follows from 
different kind of arguments. Indeed, while profile equations for genuinely-traveling waves are 
differential functional equations, profile equations for standing waves are algebraic functional 


^In a sense analogous to [?, ?, ?, ?, ?] for dissipative continuous systems. 
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equations. This results in the fact that even when traveling and standing waves coexist the 
standing-wave limit c —)■ 0 is very singular. For this reason, our analysis focused on on genuinely- 
traveling waves do not extend in a straightforward way neither to breathers nor to steady 
spatially periodic solutions. 

Our goals are, first, to derive — on a formal basis — a system of partial differential equations 
that is expected to describe the time evolution of slowly modulated wavetrains of lattice dy¬ 
namical systems, second, to elucidate — on mathematical grounds — what are, at the spectral 
level, the connections between the linearization about a constant state of the formally-derived 
modulation systems and the linearization about a periodic wave of our initial lattice system of 
equations. Our formal derivation is in many ways very classical and it is well-known as one of 
the possibilities to recover a continuous description from a discrete model. The reader is referred 
to [?, ?, ?, ?, ?, ?] for various perspectives on this point of view. Our main motivation here 
is to provide a rigorous spectral validation, that should be compared with some contributions 
of [?, ?, ?, ?] on partial differential equations, essentially as [?, ?, ?, ?] may be thought as a 
discrete counterpart of [?, ?]. 

As a direct consequence of our results stems the proof that characteristic speeds of averaged 
systems do provide group velocities for the original one. Mark that even if we illustrate our 
strategy mostly on systems obtained from semi-discretization of continuous models we provide a 
proof at the discretized level so that, in these particular cases, group velocities implicitly depend 
on the mesh size and of the particular choice of discretization. Given the fundamental role of 
a precise knowledge of group velocities in a quantitative description of the dynamics, and since 
accurate reproduction of dispersion relations obtained by linearizing about constant states has 
served for a long time now to discriminate between numerical schemes with, from other respects, 
similar performances, we do not exclude that in a near future comparison between discrete group 
velocities and continuous ones could serve similar purposes. 

Our contribution may seem relatively modest in that the proved connection is purely at the 
spectral level. Yet on the other hand our strategy is very robust as we emphasize by treating 
here in a single place various classes systems that were considered at the continuous level in 
separate contributions. Moreover, as our proofs provide relations for eigenmodes — including 
eigenvectors and not restricting to eigenvalues —, we expect that it could be directly useful for 
a dynamical validation of averaged systems in the large-time limit, as, in continuous cases, [?, 
Section 4.2] for [?, ?], [?] for [?] or [?] for [?]. 

From our point of view, our main restriction is that in the spectral validation we only consider 
periodic waves whose period is a multiple of the generator of the lattice. This stringent restriction 
does not fit well in the spirit of slowly modulated wavetrains that involve a continuum of such 
periodic waves but it enables us to use techniques originating in the study of differential operators 
with periodic coefficients — the classical Floquet theory — rather than those for quasi-periodic 
coefficients. For a glimpse at the fascinating but technically daunting quasi-periodic theory the 
reader is referred to [?, ?, ?, ?] where is analyzed the spectrum of operators similar to those 
obtained by linearizing discrete equations about a general spatially-periodic stationary solution. 

From a technical point of view, in our analysis of lattice systems, departures from the con¬ 
tinuous strategy originate in the loss of a genuine space-translation invariance of the original 
system — see Remark 3.1 though. To balance this, we make an extensive use of the remain¬ 
ing time-translation invariance and of the space-translation invariance of the profile equation. 
In particular, instead of moving to a mobile frame in which the wave is at rest and analyzing 
the spectrum of the generator of the linearized dynamics, we study directly Floquet multipliers 
of the linearized time evolution, that is, we study the spectrum of the map that encodes the 
evolution over a time 1/cj according to the linearized dynamics. Likewise, when dealing with 
Hamiltonian systems, we replace in our arguments the missing momentum conservation law with 
local conservation of energy. On the other hand, we may still use a normal, Floquet exponent 
by Floquet exponent, Bloch-wave decomposition by relying upon the discrete Bloch transform. 
Then to examine the spectrum of resulting discrete Bloch symbols we perform here a direct 
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spectral perturbation analysis ”a la Kato” [?] as in [?, ?] for continuous systems. Alternatively 
we could also have introduced a suitable Evans function, as done in [?] for continuous dynamics 
and in [?, ?, ?] for other classes of coherent structures in lattice systems. For a detailed dis¬ 
cussion of respective advantages of proofs by direct spectral perturbation or by Evans function 
expansions see [?]. 

The rest of the present paper is organized as follows. In the second section, examining 
formal expansions, we derive modulated equations. In the third section we introduce tools from 
functional analysis — including the discrete Bloch transform — that we shall subsequently use 
to analyze the spectrum of the linear one-period evolution operator. In Section 4 we state and 
prove our main results concerning the spectral validation of the previously obtained averaged 
systems. We close this paper with further comments and open questions. In Sections 2 & 4, we 
analyze three classes of systems in increasing order of algebraic complexity. First, we consider 
a lattice system of reaction-diffusion type, in the sense that the system does not support any 
local conservation law. In this case, the averaged evolution obeys a scalar equation and both 
derivation and spectral validation are considerably simpler so that we also provide high-order 
versions of those. It is actually those higher-order versions that are needed in the dynamical 
large-time validation of the slow modulation scenario. After that, we examine a general system 
that includes built-in conservation laws. The averaged dynamics is now prescribed by a genuine 
system and the spectral validation includes spectral perturbation of a Jordan block, hence 
requires a preliminary desingularization. At last, we consider a lattice system of Hamiltonian 
type that comes with an ’’extra” conservation law encoding conservation of energy. Again, all 
the corresponding theorems, Theorem 4.1, Theorem 4.2 and Theorem 4.3, are stated and proved 
in Section 4. 


Notation. Elements U G (0“^)^ of (0“^)^ are implicitly assumed to have coordinates Uj = 
U{j) G C^, i G Z. Likewise for U G Uj{t) = U{j,t), for {j,t) G Z x R, and 

U{t) = G for t G R. All along T denotes the left shift operator on (0“^)^, 

{TU){j) = U{j + 1), foTU€{&f. 

The set N is the set of nonnegative integers, and, for (M, A^) G Z, |M, N} = [M, N] n Z. The 
resolvent set and the spectrum set of an operator are denoted respectively by /)(•) and a{-). 

2. Formal derivation of averaged equations 

In this section we derive, from formal considerations, averaged modulation systems for three 
classes of nonlinear lattice systems, starting with the simplest case of lattice systems of reaction- 
diffusion type. To do so, we will follow the two-timing method introduced by Whitham [?]. Only 
afterwards, in Section 4 , shall we prove that the original dynamical evolution is indeed related 
to these modulation equations. 

2.1. Reaction-diffusion case. We consider a discretization of the scalar reaction-diffusion 
equation for u : R x R”^ —> R'^ 

(2.1) dtu = Au + f{u) 

where / : R'^ ^ R'’* is a given (in general nonlinear) smooth function. Using a centered difference 
scheme to discretize the spatial Laplacian operator, we obtain the following system of coupled 
ordinary differential equations 

( 2 . 2 ) ±u,{t) = fi[U,+i{t)-2Ujit) + Uj.i{t)] + f{Uj{t)), jGZ, t> 0 , 

where /x G R is a given fixed constant. To consider it in the abstract form of dynamical lattice 
systems, we introduce F the nonlinear operator on [/ G (C^)^ associated with / and given by 
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F{U) = {f{Uj))j^z- Then, we write system (2.2) for [/(•) = {Uj{-))j£z as 
d 


(2.3) 


dt 


U{t) = ^ (T - 2Id + T"^) U{t) + F{Uit))), t>0. 


A profile u generates through U{t) = {u{kj + wt))jgz, a periodic traveling wave to (2.3) with 
spatial wavelength k ^ 0 and time frequency oj if and only if 

(2.4) cju'iC) = fi[u{C + k)-2u{C)+u{C - k)] +f{u{C)), 
and u is one-periodic 

(2.5) u(C + 1) = ^x(C), CgR- 

To obtain (even formal) pieces of information about a given periodic wave, we need a pre¬ 
cise knowledge of neighboring periodic waves. With this respect we make the following ’’non 
degeneracy” assumption: 

Periodic traveling waves of (2.3) — identified when coinciding up to time 


(A) 


translation — form a smooth manifold of dimension one, regularly parametrized by k. 


For latter use, we denote by u^{-),io{k) such a parametrization by k. We refer the reader to 
references given in the introduction for proofs of existence results. We also stress that our as¬ 
sumption is consistent with both the continuous limit leading to (2.1) and similar considerations 
in [?, ?, ?, ?, ?, ?]. 

Our formal ansatz is of two-scale type, one scale slow — but otherwise arbitrary —, involving 
macroscopic variables, another one fast but oscillatory, hence described through an oscillation 
phase. Explicitly we want to gain some insight about the asymptotic behavior — when the slow 
frequency e goes to zero — of families of solution expanding as 





+ 0(e^), 


( 2 . 6 ) 

with phases 

(2.7) = (l)o{X,T) + e(l),{X,T) + e^cl)2{X,T) + 0{e^), 

and profiles one-periodic in the third variable 0 

u(^)(A,r,0) = uo{x,T,e) + €ui{x,T,e) + oi€‘^). 
This yields the following relation 

, a.,,-) + (e,,(..(x.T)) (v.r. 


= h 


„(.) (x +r. AlWtiA)) _ (X. T, AIAA ) + (x -T. 


+ A,T, 


4>^\x,t) 


+ 0{e^). 


Usual identification of powers of e and replacement of the phase 4>o{X,T)/e + (j)i{X,T) with 0 
then provide at order 

{dTMX,T)) deuo{X,T,Q) 

fi[uo{X,T,e + dxMX, T)) - 2uo (A, T, 0 ) + «o (A, T, 0 - dxM^, T)) ] 

+ f{uo{X,T,G)), 
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which is, in the 0 variable, the functional profile equation with wavenumber and time frequency 

( 2 . 8 ) ko{X,T)-.= dxMX,T) and iOo{X,T) := dTMX,T). 

Hence, under a suitable normalization of parametrization, without loss of generality 

(2.9) uo{X,T,e) = u*^o(^’^)(0) 
where evolution of the local wavenumber obeys 

( 2 . 10 ) drko — dx{^{ko)) = 0 , 

stemming from Schwarz’ identity 

As a result of our heuristics, we obtain the conjecture that slow/oscillatory solutions evolve 
at main order according to a slow modulation scenario 

Uf\t) = ( ^o{e3,et) 

that is, locally, at scale 1 , the solutions look like one member of the periodic traveling wave family, 
but with local parameters — phase shift (pi and wavenumber fco — evolving on a slow scale 1 /e. 
In this description local phase and local spatial wavenumber are tied by fco = dx(po and the 
evolution of the local wavenumber fco obeys at main order (2.10). In particular, unsurprisingly, 
in this simple context, linear group velocity around a given wave U_ generated by a profile u- is 
expected to be u}'{k). 

In this simple case where the family of periodic waves is one-dimensional and hence modulation 
equation ( 2 . 10 ) is scalar, algebra required to raise the order of description remains relatively 
simple so that we explain it now. It requires a more precise account of assumption (A). In 
assuming a regular parametrization by k, we include that for any wave under consideration — 
u with wavenumber k and frequency a; — the operator L acting on Lpg^([ 0 , 1 ]) with domain 
Hpg^([ 0 , 1 ]) through 

(2.11) Lu = —uiu' + iJ.[u{- + k) — 2u + u{- — k)]+df{u{-))u 

has one-dimensional kernel. By time translation invariance, this kernel is thus Cu^ Moreover 
since, as a (relatively) compact perturbation of — w( )/ L is Fredholm of index 0 — when 
considered as usual as a bounded operator on its domain endowed with its graph norm —, we 
conclude the existence of a G H^g^([ 0 , 1 ]) spanning the kernel of the adjoint C* of C and such 
that u') = 1 (in the Lpg^([0,1]) sense). In particular by differentiating the profile (2.4) and 
taking a scalar product with one receives 

( 2 . 12 ) 

Moreover, allows us to give a precise account of what is the appropriate normalization of 
parametrization that we have alluded to above. Indeed we explicitly require 

(N) {yt,dkU^) = 0 

in order to enforce the distinct respective roles of wavenumber and phase shift. The above 
condition may be achieved by shifting the original in a suitable ^-dependent way. 

After these preliminaries we may go back to the identification of terms of order e^. By 
dropping all {X, T) dependencies, we receive 

-Cui + dTpideuo - dx(pi fJ- ['9euo(- + ^o) - deuoP - /cq)] 

= -OrkodkiJ'^ + fidxko [dku^°{- + ko) - - ko)] 

+-lidxko [50Uo(- -k ko) -k demi- - ko)] 


6 



where C is here an (X, T)-parametrized version of original C, associated with uq and acting in 
the 0-variable. Averaging against yields 

8x4’! - oj'{ko) ki = d{ko) dxko ■ 

where, for concision’s sake, we have introduced a diffusion coefficient 

(2.13) d{k) = n [dku^{- + k) - dku^i- - k)] + + k) + - /c)]). 

As a refined conclusion, we obtain that slow/oscillatory solutions evolve at main order ac¬ 
cording to 

Uf{t) = ^ 

with local phase and local spatial wavenumber tied hy k = 8x4 and the evolution of the local 
wavenumber k obeys at second order 

(2.14) 8Tk — 8x{oj{k)) = e8x{d{k)8xk). 

Mark that in this second-order description, there is no more undetermined local phase shift 
and that, when present, large-time diffusive decay or anti-diffusive growth may be captured by 
equation (2.14). 

2.2. General mixed case. The main point about discrete reaction-diffusion systems consid¬ 
ered in the previous subsection is that they do not support any built-in conservation law. As¬ 
sumption A then requires that indeed no hidden conservation is present. This leads to a scalar 
modulation behavior encoded by (2.10) or (2.14). We now relax this stringent assumption to 
consider a system mixing conservative equations with non conservative ones. 

At the continuous level such structure is ubiquitous. For instance it naturally emerges from 
the modeling of the evolution of isentropic compressible flows undergoing external forces by the 
Navier-Stokes or Euler systems. Typically, in the former case one receives for the evolution of 
mass density p and velocity field u a system in the form 

8tp + 8^{pu) = 0 

8t{pu) + 8xipu‘^ + Pip)) = 8xiiyip)8xu) + gip,u), 

where P, v and g provide respectively the pressure law, Lame viscosity coefficient and external 
forces. 

Here we consider for U = (R, IT)^ — valued in (R'^)^ = (R'^i)^ x (R'’*^)^ — systems of the 
following general form 

^Rit) + B,ifrm),Wit))) = 0 

(2.16) 

^Wit) + B2iURit),Wit))) = D3(R(R(f),lT(f))D4lT(f)) + giRit),Wit)), 

where / = (frjw)^ : ^ g : ^ R'^^ ^nd B : R^ ^ CiR^^) are identified with 

their component-wise counterparts actions on sequences, of respective type (R'^)^ ^ (R'’*)^, 
(R'’*)^ —> (R'’*2)Z (R'^)^ —> CiiR^'^)'^), andDj, j = 1,2, 3,4, are constant-coefficient discrete 

differential operators, Di being conservative in the sense that the kernel of (the adjoint of 
D 1 on 4^i7j]R^^)) contains constant sequences. From now on, for definiteness’ sake, we will 
restrict our attention to the case where for some p > 0 

(2.17) Di = D4 = r,(T-Id), D 2 = Dg = 7 ?(Id-T-^). 

As in the reaction-diffusion case, we make a ’’non degeneracy” assumption on the set of 
periodic traveling waves: 

Periodic traveling waves of (2.16) — identified when coinciding up to time 
(A) translation — form a smooth manifold of dimension di -|- 1, regularly parametrized 
by (/c,M), their wavenumber and the average-values of their dith first components. 
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Explicitly M = (Mi, • • • , with Mj = ej ■ u with ej the jth element of the canonical 
basis of For latter use, we denote by a;(A:, M) a parametrization provided by the 

previous assumption. Again we refer the reader to references in the introduction for proofs of 
existence results. 

Again, a formal ansatz 


(2.18) 





+ 0(6^), 


with phases 


(2.19) <P^^\X,T) = 0o(^,T) + e0i(A,T) + e2 02(A,T) + O(e3), 

and profiles one-periodic in the third variable 0 

u^^\X,T,e) = no(X,r,0) + eni(A,r,0) + C)(e2) 


leads, under a suitable normalization of parametrization, to 

Uo{X,T,e) = ^{fco,Mo)(X,T) (0^ 

with 

ko{X,T) := dxMX,T ), uJo{X,T) := dTMX,T) 

and 

Mo{X,T):= [\ld, Od,xd,)no{X,T,e)de. 

Jo 

But now the law of conservation of waves 


drko — dx{oj{ko,M.o)) = 0 

fails to describe completely the time-evolution. 

To proceed, one collects terms of power e^, replacing 4>q{X,T)/ e + 4>i{X,T) with 0, and, with 
usual implicit notation for coordinates splitting u = {r,w)'^, receives 

-£ui + driuo) +'nTko{df{uo) ■ dxuo) = 90 ('" ) + (^.*^.^ > 


where C is the operator associated with the linearized profile equation and acts on functions 
of 0 through Tfc(,(r;) = v{- + ko). Motivated by the fact that constant functions with values 
Sj, j = 1,... ,di, lie in the kernel of £*, we average over 0 G (0,1) the scalar product of the 
previous system with those ej and obtain 

9tMo -I- dxF{ko, Mq) = 0 


where 

F(A:,M) = ri [' Mu^’^iQ)) dQ . 
Jo 

Altogether we have derived the following averaged system 


( 2 . 20 ) 


dxk — dxu}{k, M) = 0 
dTM + dxF{k,M) = 0 


for the evolution of local parameters involved in a slow modulation description of the dynamics. 



2.3. Hamiltonian case. Up to now we have implicitly assumed that our systems do not contain 
any ’’hidden” conservation law. In doing so we were motivated by our will to restrict to non¬ 
degenerate cases. However for special algebraic classes of systems ’hidden’ conservation laws are 
indeed generic and we need to accommodate them. As a typical example we analyze systems that 
are (at least formally) of Hamiltonian type. Incidentally, although we do not follow this path 
here, we stress that the presence of an extra structure also offers alternative ways of deriving 
the same modulation equations; see [?]. 

To emphasize analogy with results on continuous dynamical systems in [?] we focus on a class 
of Hamiltonian lattice dynamical systems that include discrete counterparts to the Korteweg- 
de Vries equation and the Euler-Korteweg system. Explicitly we fix some Hamiltonian % : 
(R‘’')2 ^ R and some skew-symmetric constant-coefficient discrete differential operator J in a 
conservative form, say for dehniteness’ sake J = D B with B a symmetric matrix of size d and 
D = 2 (x — . Moreover we choose some other constant-coefficient discrete differential 

operator D (not necessarily skew-symmetric), say D = 77 (T — Id), for some 77 > 0. 

Then we consider the following lattice dynamical system 

( 2 . 21 ) ^U{t) = idUpit)], 

where 5 denotes a discrete Euler operator providing variational derivatives for functionals stem¬ 
ming from local functions of ([/, DU), 

(2.22) 5H[U\ = Vt/?^(U,DU) + D*V5^?^(U,DU). 

By dehnition the previous system comes with a local conservation law for the discrete Hamil¬ 
tonian 

\T-\5n[U])B5'H[U\+T-\V^jj'H{U,Y)U)) ■ 3 5n[U] . 

Notice that to derive (2.23) from (2.21) we have used that 

A DR-(DM) H = D((T^M)H) 

for any A, H G 

As we have already done in previous subsections, to the shift operator T acting on sequences 
we may naturally associate a family of shifts on functions Tfc(7;) = v{- -|- k). Eor concision’s 
sake, from now on we will likewise denote by {P{T))k the operator P(Tfc). In particular, we 
will freely use notation D^, D^, 5k-■■ 

Now, our ” non-degeneracy” assumption on the set of periodic traveling waves takes the fol¬ 
lowing form 

Periodic traveling waves of (2.21) — identified when coinciding up to time translation 
(A) — form a smooth manifold of dimension d + 2, regularly parametrized by {k, M, E), 

their wavenumber, the average-values of their components and their Hamiltonian. 

More explicitely 

E = [ 

Jo 

With notational convention similar to the one of previous subsections, a two-scale ansatz 
leads, under a suitable normalization of parametrization, to 

UoiX,T,e) = uikoMo,Eo){X,T) 

with 

koiX,T) ■.= dxMX,T), u;o{X,T) 

Mo{X,T) := [ uo(X,r,0)d0, 

Jo 
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(2.23) ^[iU(U,DU)] = D 



and ^ 

Eo{X,T) := / H{uo{X,T,Q),T>f,^^x^T)MX,T,Q))dQ. 

Jo 

Moreover we still have a law of conservation of waves 

drkQ — 5x(<^(fco, Mo)) = 0 
and a conservation law for the averaged-values 

9tMo = BdxF{ko,M.o,Eo), 

where 

Fik,M,E) = V [" Sk{n)[u^^^^^]{C)dC = V /'vc/7^(ii'^’^’®(C),Dfcn"’^’^(C))dC. 

Jo Jo 

The latter fact is easily deduced from the fact that by collecting terms of power stemming 
from the insertion of our ansatz in our system we receive an equation in the form 

-Cui + driuo) - +T.k,)Bdx{Sko{n)[uo]) = 5©(-• •) + Dfc„(. • •) • 

The derivation of an equation for the time evolution of Eq requires more algebra. To write 
computations in a compact way we introduce some more pieces of notation. First we set := 
e“'^' then we denote and the operators involved in the expansion 

= C + (ieA;)£« + (i^A:)2/:(2) + 0(|CP). 

With these preliminaries the equation obtained from the collection of coefficients of is explicitly 
written as 

—Cui + ( 9 t ( wo ) + {dT4>i + c{ko,Mo,Eo)dx4>i) Oquo 

= dxcpiC^^^deuo + E^'^'^deuo + - c{ko, Mq, Eq)) dxuo . 

Since 5fco[^](^o) belongs to the kernel of £*, we shall average over 0 € (0,1) the scalar product 
of the previous system with SkQ[H]{uo). Our claim is that we receive 

OtEo - dx{S{ko,Mo,Eo)) = 0 

where 

Sik,M,E) = ?7 /’'iT_fc(4(^)[n('^’'^’^)])B4(W)[ll^’'^’®)] 

Jo 

+ r? [' ■ JkSkin)^^^^ 

Jo 

Indeed it does follow from direct computations. One one hand, when a G {C, k, Mi, • • • , M^, E}, 
one computes that 

(2.24) dk[n]{u) • dau = da{n{u,Uku)) - vT-k{x^^n{u,Bku)) • d^u dak + (• • •) 

where (• • •) denotes mean-free terms and dependence of u on parameters has been omitted. On 
the other hand, for a G {C, k, Mi, • • • , M^, E}, we also have that 

(2.25) 

4['H](w) • - c)daU + 

= da(p\T_kmn)M)B5k{nM + r]T_k{Vf,u'^{u,Bku))Mk8k{UM) 

- 'nT_k{^^jjU{u,Bky^) ■ di^udaUJ + (■■■) 

with the same convention. This yields our claim by integration using that drko = dx(uj(ko,Mo, Eq)). 
Alternatively one could have derived this last equation by directly expanding and averaging 
equation (2.23). 
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Altogether we have derived the following averaged equations for the slow modulation of local 
parameters 


( 2 . 26 ) 


' drk = dxio{k,M,E) 

< drM = BdxF{k,M,E) 
, OtE = dxS{k,M,E). 


3. Analytic framework 

We collect in this section some preliminaries from functional analysis that will allow us to 
formulate in which sense the formally derived averaged equations provide some valuable pieces 
of information at the spectral level. 

3.1. Integral transform. To analyze linearized lattice systems with periodic coefficients we 
shall make use of an adapted integral transform, usually called discrete Bloch transform. Though 
probably less known that its continuous counterpart it is a natural functional tool that have 
already received various applications in related contexts, see for instance [?, ?]. 

3.1.1. Discrete Bloch Transform. Let N G N* be given. We define the following transform 

(3.1) B: ^ l2(|0,7V-11x , / ^/ 

where for j £ |0, A^ — 1] 

(3.2) f{j,-) := lim ^ ' fkx+j, 

ko^oo 

A:€[—fcoTol 

the limit being taken in the ([~;^) ;^])-sense. Up to a multiplicative constant this is a total 
isometry whose reciprocal is given as 

(3.3) B~^ : (|0, iV - 1] X [-f, f ^ i\Z- C^), f f 

where for {k,j) G Z x |0, — 1] 

(3.4) hM+j = L 

From now on we’ll constantly identify functions on Zf{NZ) with functions on |0, A^ — 1] 
through usual restriction/periodic extension processes. To this purpose we stress that for¬ 
mula (3.2) already provides the needed periodic extension. With this identification, inverse 
formula (3.4) simply reads for any J G Z 

(3.5) /, = — / ^e‘^«/(i,C)dC. 

3.1.2. Discrete Bloch symbols. To any element ^ i-7> of ; C{i‘^{Z/{NZ)-,C'^))) 

one may associate a bounded operator A G C{£^{Z; C^)) such that, for any / G ^^(Z; C'^), 

iAfn;0 = 

Explicitly 

(3.6) ( Af ), = ^ U (.45/(.,0),d{. 

Moreover up to a multiplicative constant the map 

L-([-f,f] ■C{f{Z/{NZ)-C‘^)))^£if{Z-C<^)) 
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is a partial isometry. From this follows the classical identity between the spectrum of the 
synthetized operator A and the essential range of i—> a{A^). In particular, in the case where 
Ag belongs to {NZ)-,C<^))), this yields 

= U 

so that the spectrum of A may be analyzed by looking at the family of spectra of the finite¬ 
dimensional symbols A^. 

Reciprocally for a large class of G C{£‘^{Zi; C^)) one may define corresponding symbols 
^ I— )■ formally given by 

A^ = e“‘'^ yl e^'^ 

that is, for / G f{Zj{NZy,C'^) 

(A (e‘-V)), 

with an extended sense of A. To make things more concrete, we focus now on the case where A 
is an A^-periodic discrete differential operator in the sense that it is given as 

A = a(-,T) with a{-,X) = 

where T is again the left shift operator and a ; Z —>■ C{C^) is Ai-periodic. The Bloch symbols 
of A are then given by 

= a(-,e'«T) 

where T is now the left shift operator acting on and a is considered as a function 

on ZjiXZ). 


3.2. Linear evolution. By linearizing original equations we receive equations in the form 

iF(t) = A(t)V{t) 

where yl(-) is a 1/w-periodic family of bounded operators generating 5(-, •) an evolution system 
on l‘^{Z-,C^) (in the sense defined for instance in [?, Chapter 5]). Explicitly we have 

(Vs G R, S{s,s) = Id) and (V(s, r, t) G R^, S{t,s) = S{t,r)S{r, s)) 

and 

V(s,t)GR^, {dtS){t,s) — A{t)S{t,s) = 0 and {dsS){t,s) + S{t,s)A{s) = 0. 
Mark that periodicity implies that linear stability in the sense that 


sup ||5(t,0)|| < oo 
teR+ 


is equivalent to 


sup ||5(l/w,0)* 

fceN 


< oo. 


Likewise linear exponential instability in the sense that 

37>0, inf e“'*'* ||S'(t,0)11 > 0 

t€R.+ 


is equivalent to 


3k > 1, 


inf 

fceN 


k-"||5(1M0)" 


> 0 , 


that is to p(5(l/w, 0)) > 1. 


12 



Remark 3.1. Here we use time-translation invariance and periodicity in time as a substitute to 
space-translation invariance and stationarity in a moving frame that play a crucial role in similar 
analyses of continuous systems. However as pointed out for instance in [?, ?, ?, ?] the loss of an 
exact invariance in space does not preclude arguments based on a pseudo-invariance. Notably, 
mark that for all t G R'*' 

A{t + l/c) =T-M(t)T. 

This implies that for all {t, s) G R 

S{t + 1/c, s + 1/c) = s)T . 

As a consequence, linear stability is also equivalent to 

sup ||(T5(l/c, 0))^|| < oo . 

fceN 

The fact that such an alternative exist is obviously crucial when dealing with fronts, shocks or 
solitary waves rather than with periodic waves. 

From now on, in order to make the most of the discrete Bloch transform we focus on cases 
where, for some N G N*, A(-) is valued in N-periodic discrete differential operators that is we 
restrict to cases where linearization is about a periodic traveling wave with an integer period, N. 
Then A^(-) the value at ^ of the corresponding family of Bloch symbols generates the evolution 
system given by the value at S^{-, •), of the family of Bloch symbols associated with S{-, •). 
In our context some continuity in the ^-variable is available and linear stability reads 

sup sup ||5^(l/w, 0)^11 < oo , 

N'Nl 

while linear exponential instability may be written as 

3Cg , p(5^(l/a;,0))>l. 

This implies that a necessary condition for linear stability is that for any ^ G 

C R(0,1). 

The fact that each periodic wave is embedded in a family of neighboring waves implies that 1 
always belongs to the spectrum of 5o(l/w, 0). It is then crucial in the analysis of linear stability 
to determine how the eigenvalue 1 evolves when ^ is varied. Our goal is to prove that this piece 
of information is precisely given by the averaged equations. 


3.3. Algebraic interplay. To analyze the evolution of the spectrum of 5'^(l/a;,0) when ^ is 
small, we first stress the following classical identity, for any ^ G ;^] 

(3.7) V(s,t)GR^, S^{t,s) = So{t,s) + So{t,r){A^{r) - Ao{r)) S^{r,s)dr. 

In particular, expanding Ag(-) as 

A^{-) Ao(-) + iC A«(.) + (iO^ + 0(|^|3), 

uniformly in time, for some we receive locally uniformly in time 


with 


(3.8) 


5^(t,s) 5o(t,s) +ie5«(t,s) + (iO' 5(2)(t,s) + O(ieP), 

= J 5o(t, r) A^^^(r) 5o(r, s) dr 

= f 5o(t, r) A^^^(r) 5o(r, s) dr -|- f So{t,r) A^^\r) s) dr . 

J S J S 
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We have already encountered a similar expansion for := Operators acting 

on functions of period 1 are actually continuous Bloch symbols of the operator obtained by 
extending C to functions defined on the full line. Since we won’t need this relation here, we 
mostly use the fact that by expanding 

= C + (iCA:)/:« + (i^A:)2/:(2) + 0(|^|3). 

we obtain operators that are practical shorthands to denote relations obtained by differen¬ 
tiating profile equation with respect to parameters. For instance, when a denotes a parameter 
which is not k (but possibly (, to account for phase shifts), one receives 

(3.9) CdaU = daOjdc^u 

while 


(3.10) 


CdkU = —kdkcdc^u — C^^^d(^u. 


Proofs of our main results strongly rely on the interplay between and that allows to 
use the previous relations. Explicitly to any function u : R —of period 1, we associate 
W : R ^ (C'^)^ by 

E^(j,t) = v{k{j-ct)). 

Observe that when v and w are one-periodic 


^'^‘"(E-(s),W(s))ds = ^{w,v). 

Note also that when v is one-periodic is 1/w-periodic in time and valued in Wperiodic 

sequences. In this case, 

- A^{t)V^{t) = 

^W(t) + A^{t)*V^{t) = E^«^(t) - iewW(t). 

In particular 

(3.12) A((F-,W))(t) = - (E-(t),E^«"(t)), 

which turns to be crucial to transport duality relations. Obviously by expanding in ^ the 
foregoing relations we obtain corresponding relations between operators A^^^ and operators 
for instance 


(3.13) 


Aw(t) - Aoit)V-{t) 

- A^^\t)V^{t) 

- A(2)(t)E«(t) 

Aw(t) + Aoityv^t) 


Mark that from (3.9) and (3.11) 

V(t, s) G R2 , = 




So{t,s) V^c^s). 


In particular, E^^-(O) G ker(So(l/ui,0) — Id). More generally (3.9) and (3.11) provide as many 
members of the generalized kernel of S'o(l/w, 0) — Id as the expected dimension of the family of 
periodic traveling waves (counted up to translation). Regardless of the version of the equations 
that we are considering, assumption A may be reformulated as the fact that these elements form 
a basis of this generalized kernel. Moreover (3.10), (3.11) and (3.8) yield 


(3.14) V(t,s) G R^ 


So(t,s)V^'‘^(s) - (t-s)dkUjV^c^(s) + i5«(t,s)E^c-(s). 
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4. Spectral validation of averaged equations 


4.1. Reaction-diffusion case. Firstly we analyze the linearization of (2.3) about U_{') the 
periodic traveling wave associated to n-, with k = 1/N. The linearized evolution obeys 
(Vt, = ^(t) y(t)) where 

T(.) = M(T-2Id + T-i)) +dfm-)) 

with usual identification of functions with their local action on sequences, that is, here, 

v(t,j), [dfmtWh = df{U^{t))V,. 

Explicitly here 

= and {T + T-^). 

Now our main goal is to prove the following result that validates (2.10) and (2.14) at the 
spectral level. 

Theorem 4.1. Assuming that condition A holds in a neighborhood of there exist positive cq 
and .^0 G (0)7r/N) and an analytic curve A ; [—C such that 

V^€Ko,^o], ct(55(iMo)) n R( 1 , 6 o) = {A(0}. 

Moreover, for f G [—Co)?o]; eigenvalue X{f) of S^{l/uj-,0) is simple and 

(4.1) A(0 0(|^P), 

where d{k) is given by (2.13) under normalization N. 

The above theorem has many direct implications. For instance, we readily see that d{k) > 0 
yields a linear exponential instability caused by side-band perturbations, that is, by pertur¬ 
bations with arbitrary small — but non zero — spatial Floquet exponents. It also proves as 
expected that the single linear group velocity is indeed dkU)-. 

Proof. Assumption A implies that 1 is a simple eigenvalue of 5o(l/w, 0) so that the only point 
that does not follow from classical results on regular perturbations of simple eigenvalues is 
expansion (4.1). To identify coefficients we expand 

A(6 1 + iC A« + {iff A(2) + 0(|e|3), 

and a corresponding eigenvector which also depends analytically on f, as 

f9cm( 0) + <^(1) + (i^)2 (^(2) + 0(|^|3). 

Setting p = 1/uj- for writing convenience, we now expand relation S^{p,0)ip^ = X{f)^p^. 

At first order in i^ we receive 

(4.2) (5o(p,0)-Id)(^W = aW ■F^c^(O)-5W(p,0)F^c^(0). 

Hence 

{Soip, 0) - Id) ((^(^) - kV^>^^0)) = (a(i) - tpdkuj^'j F^«^(0). 

Since 1 is a simple eigenvalue of Sq{p,0), F^^-(O) does not belong to the image of So{p,0) — Id. 
From this we infer the desired 

A(i) = l^dkio^ and - kV^’^^O) G kefSaip, 0) - W) = CF^f^(O). 

By expanding now at the second order in i^ we obtain 

(4.3) (5o(p,0)-Id)^(2) ^ a( 2) -^^^1(0)-5W(p,0)(/jW - 5(2 )(p^0)F^c^(0). 

To proceed we also need to involve defined in Subsection 2.1, that is spans ker£* and 
iiLk^j d^^vA) = 1 (in the Lpg^([0,1]) sense). From (3.11) we deduce that 

V(t,s), So{s,trV-^{t) = V^{s), 
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-,ad 

thus that (0) G ker(5o(p, 0)* — Id) and from (3.12) that 

Vt, = - [-{V^t(s),V^<^s))ds = l{uf,dcv^) = i 

PJo H H 

Since —^dkUJ-V^<-{0) + 0)F^«-(0) belongs to the range of So{p, 0) — Id this yields 

^A(2) = (V-l\o),-^dku;^V^^^(0) + kS(^Hp,0)V^>^^(0) + S(^^(p,0)V^<^(0)}. 

Now with (3.8) we compute 

ad 

(y-fc (0) , fe5(i)(p,o)y^fcii(o) + 5(2)(p,o)y^cii(o)) 

(F-i‘'(s),A(i)(s)[fe5o(s,0)y^'=2^(0) + 5(^)(s,0)y^ci‘(0)] + (s)V^<^(s)} ds 

(V-^\s),A(^^(s)lksdku;^V^‘:^(s) + kV^'^^(s)] + A^^^(s)V^(^(s)} ds 

Now, since (£(^^ — c)d(^u = dkOjd(^u — Cd^u^ by using (3.13) and normalization N and integrating 
by part, we observe that 

\ky{dkOJ-? + k\dkUJ^ {V^t (0), y^'=2^(0)) 

which yields 

= pl^ d{k) + ^ {pkdkioy 

as desired. □ 

4.2. General mixed case. Now our main goal is to unravel the role of (2.20) in the spectral 
analysis of the linearization of (2.16). 

Theorem 4.2. Assuming that condition A holds in a neighborhood ofu-’—, there exist positive 
eo and G (0, tt/N) and eontinuous curves Aq : [—^o] C, a = 1, • • • , di + 1, sueh that 

V^gKc^o], c 7(5^(1/^-’-,0)) ni?(l,eo) = {Ai(0,--- ,Adi+i(6}- 
Moreover, for G [—^o] o,nd a = 1, • • • , di + 1 

(4.4) A„(0 ^=° + O(ieP), 

where ai,-- - ,0^1+1 are the eharacteristie speeds of the Whitham system (2.20) linearized about 
parameters (k, M). 

In contrast with what happens in the reaction-diffusion case, that is averaged to a scalar equa¬ 
tion, the above first-order expansion already provides an instability criterion. If the Whitham 
system (2.20) is not weakly hyperbolic at (^, M), that is, if some of the Oa are not real, then the 
corresponding wave is linearly exponentially instable to side-band perturbations. In contrast, 
when weak hyperbolicity is met, then the above theorem proves that those Oa do give the di -|- 1 
linear group velocities, a non trivial fact. 

Proof. Let us denote by ei, • • • , both the di first members of the canonical basis of and 
the constant functions with corresponding values. As a consequence of assumption A, there 
exists uf^ belonging to the generalized kernel of C*, orthogonal to dMiH) • • • and such 

that {uf^,d(^u) = 1. Moreover, by appealing to Subsection 3.3 as in the reaction-diffusion case, 
one obtains that - ■ ■ ,1/^’^rfi-) and ‘ ,kV‘^'^i) form dual bases of 

the generalized kernels of 5o(l/a;-’—, 0) — Id and 5o(l/a;-’—, 0)* — Id. 
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r{V-t{s),k^sdk 0 j!^V^^‘'"^-^'>^<^s)) ds = y 

Jo k. 




By Kato’s perturbation method [?, pp. 99-100] we may continue these bases as dual bases 
:Qdi{0) (?o(0r'' yQdiiO) of fho generalized eigenspaces of ,0) and 

5^(l/a;-’—, 0)* associated with their spectra in 5(1, eo) as long as j^j < provided and cq 
are small enough. Then, for j^j < the spectrum of 5^(l/a;-’—, 0) — Id in 5(0, eo) is the 
spectrum of the matrix 


ric = 


QaiC), ('S'g(l/‘^-’—,0) - Id) q/si^) 


J 0<(y,^<di 


Observe that 


On — 


/O ••• 

VO 


where p = l/cj-’—. Expanding 


0^ = 00 + iMO(i) + (iM)^0(2) + C>(|^|3), 


we also obtain since 

i5(^)(p, 0) E^c-(o) = -(5o(p, 0) - Id) E^'=^^(0) +pdkCo!^ V^<^0) 


that 


= 


/* 

0 

Vo 


The upshot is that the matrix defined by 

Of := ——Sr^OcSf with 

^ iM ^ 


/ 




(iM)-' OixdiV 

Odj^xl Idixdi/ 


depends analytically on The heuristics behind this transformation is that the change of basis 
transforms a phase-like coordinate to a wavenumber-like coordinate as in the derivation of the 
Whitham system (2.20) while the division by ik^ changes eigenvalues in velocities. 

Now the point is to identify Oq. We already know that 


Oq 



(^o)a,/3 

~W^)^ 



0<a,l3<di 


and, for /3 = 1, • • • ,di, (Oo)i ,/3 = First 


(O0))o_o = ^(go(0), (5'o(p,0) - Id) 5^go(0)) + |(go(0), O)go(O)) 


= pdkLV + i|(go(0), {Sq{p, 0) - Id) {d^qo{0) - ikV^'^-{0))). 

Our claim is that the last term vanishes. Indeed since {So{p, 0) — Id)go(O) = 0 and, for j^j < 

(S’^(p,0) - ld)go(0 e Span({go(0r • • ,9di(0}) 

we conclude 

^{Soip,0)-ld)d^qo{0) + lS^^\p,0)qo{0) G Span({go(0), • • • , 9di(0)}) 

hence 

(5o(p,0) - Id) (5^go(0) - ifeE'^'=-(0)) G Span({go(0), • • • ,gdi(0)}). 

This implies 

d^QoiO) - ifeE^''-(0) G Span({go(0), • • • ,qdi (0)}) 

which in turn yields 

(5o(p,0) - Id) (d^qoiO) - ifeE^'=-(0)) G Span({gi(0), • • • ,gdj(0)}), 
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proving the cancellation. To make the remaining computations easier we pick (zi, - ■ ■ ,Zdi) such 
that 

l<Q:<dl 

and replace (go(0> ''' ,9di(0) with {qoiO - Zaqa{0,qi{0 + ■ ■ ■ ,9di(0 + 

^Zdj^qoiO)- This ensures d^qo{0) = ikV^>^-{0). 

With this in hands we may also compute, for 1 < a,/? < di, 

(5o(p,0) - Id) 9/3(0)) + l{qa{0),S‘'^\p,0)qfi{0)) 

= ^^-(5€9a(0),9o(0)) + |(ga(0),5'(T(p^0)9/3(0)) 

= -^^f-(9a(0),5^9o(0)) + |(ga(0),5T)(p, 0)9/3(0)) 

= -pdMfsU}{qa{0),V^'^^{0)) + |(9a(0),5'(^)(p,0)9/3(0)). 

Now using that, for any {s,t), 

Soit,s)V^^-{s) = + {t - s)dMpUjV^<^t), 

we go further to obtain 

i(9«(0),5«(^,0)9/3(0)) = [-{V^-{t),A^^\t)[V^^Ht) + tdM,u;V^<m) dt 

Jo 

with 

f {V^‘^{t),Ad'){t)V^'^^-{t)) dt = p(ea, _ c)aM^;u) 

Jo 

= p [ pea- df{u{C))dMpuiC)dC 

Jo 

= -pdMpFa{k,M_) 

and, by using (£^1 _ c)d(^u = dkOjdc^u — CdkU and integrating by part, 

l~t{V^-{t),Ad){t)V^ci^{t)]) dt = -k dt 

Jo Jo 

= k j-t (W“(t), (A _ A,{t))) dt 

Jo 

= fcp(I^®“(0),F^''^(0)). 

Hence, for 1 < a,/I < di, 

= -pdMpFa{k,M)- 

Lastly we compute, for 1 < a < di, 

= ^(9^da(0),(S'o(p,0)-Id)a^9o(0)) + ^(9a(0),5'(^)(p,0)5^9o(0)) 

+ ^(^^da(0), 5 'T)(p, 0) 9o(0)) + ^(9a(0))0) 9o(0)). 

On one hand, integrating by parts, 

-p (5«da(0),5'T)(p,0)9o(0)) = {d^qa{0),So{p,t)J^^^Ht)V^<^t))dt 

= S [~{d^qaid), So{p, t)dkU} dt + A f~{d^qa{0), So{p, t) {A - Ao{t)) dt 

-Jo -Jo 

= ^ pdkOJ (d^qaio), H^«-(0)) - ^{d^qa{0), (5o(p, 0) - Id) 5^9o(0)) • 
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On the other hand, using (3.14), 

^ (9q( 0), [i5'(i)(p,0)5^go(0) +5'(2)(p,0)go(0)]) 

Jo ^ - Jo 

= p{ea,{CJ^^ - ^dku + dc_^ + pdku ^(g«(0), 5^go(0)) • 


Hence 

= p{ea,{CJ^'^ - ^dku + C^‘^'>dc_%^ = -pdkFa{k,M) 

so that as expected 



□ 


4.3. Hamiltonian case. Now we turn our attention to the role of (2.26) in the spectral analysis 
of the linearization of (2.21). Our analysis hereafter follows closely the proof of the foregoing 
subsection. The only difference stems from the presence of one extra ’’hidden” conservation law. 

Theorem 4.3. Assuming that condition A holds in a neighborhood o/u-’—, there exist positive 
eo and £ (0, tt/N) and eontinuous eurves Aq, : [—^o] C, a = 1, • • • , d + 2, such that 

V^e Ko,eo], ^(5 €(iM’-,0)) n i?(l,eo) = {Ai(0,--- ,Arf+2(0}- 
Moreover, for ^ G [—? 0 ) ^o] (ind a = 1, • • • , d + 2 

(4.5) A„(0 ^=° + O(ieP), 

where ai,--- ,ad +2 are the characteristic speeds of the Whitham system (2.26) linearized about 
parameters (k, M. E). 

Proof As in the previous subsection, perturbatively we may obtain dual bases (?o(0r '' > dd+i(C)) 
and {qoiOr" Md+iiO) of generalized eigenspaces of 5^(l/w-’—, 0) and 5 ^( 1 /cjAM>.e^ q^* 
associated with their spectra in B{1, cq) as long as |^| < provided and cq are small enough. 
Here however we start from 

{qo{ 0 ),--- ,qd+im = 

and 

{qo{0),--- ,qd+m) = {kV^^\kV^\--- ,kV^'^,kV^>^^^), 

where belongs to the generalized kernel of C*, is orthogonal to dMiH:''' ,dMalL,dEU and 
satisfies {uf^,d(^u) = 1. Then, for |,f| < ^q, the spectrum of 5^(l/a;-’—, 0) — Id in H(0, cq) is 
the spectrum of the matrix 


The cancellation discussed in the previous subsection shows that the matrix defined by 




s. = 


(iM) ^ Oixd+1 
Od+lxl Id+lxd+1 
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depends analytically on The point is again to connect with our averaged system, here (2.26). 
Manipulations as above already yield 


1 

/ dkUJ 

cImiW 

dMa^ 

Beoj \ 

II 

0 

>C1 

-BdkF 

-B^,F •• 

■ -BdM,F 

-BSeF 





^(^^^^^)<i+i,d+i ) 


where p = l/cj-’—, for a = 1, • • • ,d + l, 

(5’o(p,0) - Id) 5^go(0)) + ^(gd+i(0),5'(^)(p,0)5^go(0)) 

+ 0) go(0)) + ^(^d+i(0), 0) go(0)), 

and, for /3 = 1, • • • , d + 1, 

= ^(9^%+i(0), (5o(p,0)-Id) g/3(0)) + |(grf+i(0),5(^)(p,0)g/3(0)). 

Moreover we may enforce d^qo^O) = ikd^u. 

For /3 = 1, • • • , d + 1, setting 


0/3 


if 1 < /? < d 
E if ,5 = d + 1 


and proceeding as in the previous subsection we derive 

= I ( 559 ^+ 1 ( 0 ), go(0)) + p -c)5a^u) 

+ kda,u y^'=^(t)]) dt. 

do 

Now using (2.24) and (2.25) we also obtain 

p{5kH\u\,{C^^'^ - c)da,u) = pda,S + kda.OJ j-,V^^^t)]) dt . 

Jo 

From this stems as expected 

(n(^))d+i,/3 = [{d^qd+i{0),qo{0)) + {qd+i{0),d^qo{0))] + pda^S = pda^S. 

Lastly, on one hand, we still have 

-p{d^qd+i{0), 5'(i)(p,0)go(0)) 


J_ 


pdkUJ {d^qd+i{0),qo{0)) 


^{d^qd+i{0), {So{p, 0) - Id) 5^go(0)) 


On the other hand, 

^ (9«(0), [Y5'(^)(p,0)5^go(0) + S^‘^'>{p,0)qoiO)]) 

fP 


(F'=“(t),AW(t) F^'=2i(t)+t4w'F^^-(t) )dt + i / {V^<^{t),A^^\t)V^‘^Ht))dt 
Jo ^ - Jo 

= p{ea,{C^^^ - c)dku + C^'^^dt^u) + kdkoj [ t - vlo(t)) F^''-(t)]) dt 

Jo 

= pdkS + pdkUJ ^{qa{0),d^qo{0)) ■ 

Hence as expected (0(^))(/+i^o = pdkS. 
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This achieves the proof of 
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that yields the theorem. □ 


5. Conclusions and Remarks 

We have proved that the slow modulation ansatz captures accurately the essential spectral 
features of the space-time low^ Floquet-multiplier evolution linearized about a periodic wave 
whose period belongs to the lattice. As already mentioned this opens at least three classes of 
questions. 

In many cases one also expects that in the large-time the dynamics is effectively reduced to 
a slow modulation evolution. Our analysis provides the spectral background — including, as 
appears from our proofs, expansion of critical eigenfunctions — for a proof of such claims, either 
in the linear regime as, for continuous systems, in [?] for the Korteweg-de Vries equation, or 
in the complete nonlinear regime, as in [?] for general parabolic systems of partial differential 
equations. 

With this in mind, one expects to be able to identify on averaged systems the key features 
of the discrete dynamics. It is then natural to ask if one may design discrete systems leading 
to desired modulation properties, either with applications to the analysis of numerical schemes 
where the goal is to preserve the slow modulation asymptotics at the discrete level, or in the 
effective conception of smart materials, a fast-growing field of investigation. 

At last, even for the direct problem at the spectral level, remains the technical challenge of 
dealing with waves of general period, a seemingly quasi-periodic problem. 
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